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We report a theoretical study of the macroscopic quantum dynamics in spatially extended Joseph- 
son systems. We focus on a Josephson tunnel junction of finite length placed in an externally applied 
magnetic field. In such a system, electromagnetic waves in the junction are excited in the form of 
cavity modes manifested by Fiske resonances, which are easily observed experimentally. We show 
that in the quantum regime various characteristics of the junction as its critical current I c , width 
of the critical current distribution er, escape rate F from the superconducting state to a resistive 
one, and the time-dependent probability P(t) of the escape are influenced by the number of photons 
excited in the junction cavity. Therefore, these characteristics can be used as a tool to measure the 
quantum states of photons in the junction, e.g. quantum fluctuations, coherent and squeezed states, 
entangled Fock states, etc. 

PACS numbers: 03.67.-a, 03.75.Lm, 74.40.+k, 42.50.Dv 
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Great interest is currently attracted to experimental 
and theoretical studies of macroscopic quantum phenom- 
ena in diverse Josephson systems yj, y, |3| ■ Most of such 
systems contain just one or few lumped Josephson junc- 
tions. At low temperatures, quantum-mechanical effects 
such as macroscopic quantum tunneling, energy level 
quantization, and coherent oscillations of the Joseph- 
son phase have been observed 0, The important 
method allowing to study the macroscopic energy lev- 
els in Josephson coupled systems is the microwave spec- 
troscopy [3) S] ■ A s the frequency u) of an externally ap- 
plied microwave radiation matches the energy level sep- 
aration, one can observe resonant absorption and Rabi 
oscillations due to the population of excited levels. 

An interaction between Josephson systems and quan- 
tized electromagnetic fields opens new frontiers in re- 
search. Similar type of interaction arises between an 
atom and electromagnetic cavity modes in quantum elec- 
trodynamics (QED). Recently, such QED-inspired exper- 
iments have been performed with superconducting charge 
qubits [4j and flux qubits [a, [6| coupled to on-chip res- 
onators. In these experiments, the quantum degree of 
freedom of a lumped Josephson system couples to the 
quantized electromagnetic field of a superconducting cav- 
ity located on the same chip, as illustrated in Fig. 1(a). 
The excitation of cavity modes CMs in an external cav- 
ity leads to an appearance of an oscillating current flow- 
ing through the Josephson circuit. Such an oscillat- 
ing current excites transitions between macroscopic en- 
ergy levels of the Josephson phase. Quantum regime 
of a weak resonant interaction between the Josephson 
phase ip and CMs is described by a bilinear Hamilto- 
nian H oc ip(a + a + ), where a + (a) are the operators of 
creation (annihilation) of a particular CM. This Hamilto- 
nian corresponds to the famous Jaynes-Cummings model 
0, 0] and therefore, a quantum regime of a lumped 
Josephson circuit incorporated in an external transmis- 
sion line can be mapped to a problem of a single atom 



weakly interacting with CMs. In this regime, many fas- 
cinated phenomena as mixture of different Rabi frequen- 
cies, creation of entangled states of CMs, a single atom 
maser behavior, can be observed Q. The interaction 
via CMs can be used to couple superconducting qubits 

As we turn to spatially extended Josephson systems, 
e.g. long Josephson junctions, Josephson junction par- 
allel arrays and ladders, the spatially-dependent Joseph- 
son phase ip(x) can also display macroscopic quantum 
effects as tunneling and energy level quantization. For 
a Josephson vortex trapped in a long Josephson junc- 
tion, these phenomena have been studied theoretically 

An additional, 
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[131 ] and observed in experiments 
interesting property of spatially extended junctions is 
that they can support the propagation of electromag- 
netic waves forming CMs inside the junction itself, see 
Fig. 1(b). When the Josephson system is biased in the 
resistive state the resonant interaction of such CMs with 
ac Josephson current leads to the classical resonances in 
the current-voltage characteristics known as Fiske steps 
liL 16| . The quantum regime of such resonances at fi- 
nite voltages in systems containing few small Josephson 
junctions has been theoretically considered in Ref. [l7j ]. 
The next step naturally appearing in a study of macro- 
scopic quantum phenomena in spatially extended Joseph- 
son systems is the quantum regime of interaction between 
the Josephson phase and intrinsic CMs of the system. 

In a current-biased Josephson junction, the spatially- 
averaged Josephson phase difference ip oscillates at 
Josephson plasma frequency io p = ui p o [l — (I/I c q) 2 ~\ 
depending on the dc bias /. Here, uj p o and I c q are the 
plasma frequency at zero bias and the nominal value of 
the critical current of the junction, respectively. If the 
Josephson plasma frequency u> p matches the frequency of 
CMs, the interaction between the Josephson phase and 
CMs becomes resonant. 
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FIG. 1: (a) Schematic view of lumped Josephson junction em- 
bedded in an external cavity and its equivalent scheme, (b) 
Distributed Josephson junction subject to an external mag- 
netic field with excited internal (Fiske) CMs and its equivalent 
scheme. 



In this paper, we consider the alternative, non-resonant 
interaction between the Josephson phase tp and intrinsic 
CMs of the junction. Such an interaction can be realized 
in long Josephson junctions or their discrete versions - 
parallel Josephson arrays, often also called as Josephson 
transmission lines. We show that the switching from the 
superconducting state to a resistive one, i.e. the Joseph- 
son phase escape can be facilitated or suppressed by ex- 
citation of cavity modes. 

Moreover, we argue that measurements of the switch- 
ing current I c , the width of its distribution a, and the 
time-dependent probability P(t) of the escape allow to 
obtain information on the quantum state of intrinsic CMs 
in the junction. 

In order to quantitatively analyze the macroscopic 
quantum phenomena appearing due to the interaction of 
plasma oscillations with CMs, we consider a Josephson 
junction of a finite length L. The Josephson junction 
is characterized by the time and coordinate dependent 
Josephson phase ip(x, i). In the presence of an externally 
applied magnetic field the Josephson phase can be writ- 
ten as 



(1) 



where <J> is the external magnetic flux, and <i>o = hc/2e 
is the flux quantum [l8j]. The system is described by 
Hamiltonian: 



dc bias i — I/I c o can be changed to tune the effective 
potential U(x, (p) of the junction. 

Next, we represent the time-dependent Josephson 
phase by a sum of electromagnetic CMs and a "center 
mass" Josephson phase i.e. 



<p(x,t) = x(t) + ^2Q n (t)cos(k n x), 



(3) 



where the wave numbers of cavity modes k n = nn/L, 
with n — 1,2.... Substituting expression ([3]) in Eq. (|2|) 
and assuming that the amplitudes of cavity modes Q n (t) 
are small (so that we can neglect anharmonic interaction 
between CMs (HI), we obtain the Hamiltonian in the 
following form: 
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E 



Ho 



2 ^ 



Qn 



+ U{ X ,Qn), 



U(x,Qn) = EjoY, 

— Ej{H) cos 



x- «X , 
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where uj n = Lu p oX,jk n . The magnetic field dependent cou- 
pling energy is Ej(H) = Ej Ia j H ^ , where I C (H) is the 
magnetic field suppressed nominal value of the critical 
current, i.e. 



I c (H)=I c0 
The ^-dependent coefficients 
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determine the strength of interaction between the CMs 
and the center mass Josephson phase. 

As Ej(H) 3> fiujpQ, the switching from the supercon- 
ducting state to a resistive one occurs at the dc bias close 
to its critical value, i.e. S = [I C (H) — I]/I C (H)<^ 1, and 
the quantum-mechanical Hamiltonian ([4]) can be written 
in the following form: 



H = Hn 



where the first term 
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describes the dynamics of homogeneous (plasma) oscilla- 
tions of the Josephson phase with momentum operator 
P x . The second term is the Hamiltonian of noninteract- 
ing CMs 



where Aj and Ejq are the Josephson penetration length 
and the Josephson energy, accordingly. The normalized 
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where Pq = —ifo-g^ is the operator of generalized mo- 
mentum . Here, we have neglected a small renormal- 
ization of CMs spectrum due to the presence of plasma 
oscillations of the Josephson phase x- The last term in 
([5]) describes an interaction between the Josephson phase 
and CMs: 



Hint — — E.lX ^ a nQn 
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where a n = \j 77 cos 2 (fc ra x) cos( 2 £®£ ). In the absence of 
an external magnetic field all coefficients a n — 1/4. How- 
ever, in the presence of a magnetic field characterizing by 
magnetic flux, $ ~ m$o, the coefficient a m = —1/8, and 
other coefficients are small. 

If the size of the Josephson system is not very large 
(L < Aj), all frequencies u> n lo v q and, therefore, 
the interaction between the CMs and the plasma os- 
cillations is non-resonant. Thus we can use a time- 
averaged expression for the interaction Hamiltonian, 
H^t = EjxY,n a n < Q n >> where <_... > is the time 
average. The mean switching current I c is written as 



I c = I C (H) - Icq a„ 



<Ql > 



(9) 



The excitation of CMs results in a nonzero value of 
< Q n >, which leads to a suppression of the average 
switching current in the absence of magnetic field. How- 
ever, when the magnetic field is applied, the fluctuation- 
free critical current I C (H) is strongly suppressed and ex- 
citation of CMs results in an unusual increase of the 
mean switching current [a n values become negative in 
this case). 

Switching from the superconducting (zero-voltage) 
state to the resistive state occurs at random values of 
the external current I. This process is characterized by 
the probability P(I). If we neglect the CMs, the Joseph- 
son phase escape in the quantum regime is determined 
by macroscopic tunneling, and the width a of the distri- 
bution P(I) strongly depends on magnetic field: 
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Quantum fluctuations induced by zero-point oscillations 
of CMs lead to a saturation of a as a function of magnetic 
field at the level 
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(11) 



The typical dependence of the width of the critical cur- 
rent distribution on magnetic flux applied to the junction 
is shown by solid line in Fig. 2. 

In the following, wc present a quantitative analysis of 
the influence of quantum-mechanical properties of CMs 
on the process of escape from the superconducting state 
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FIG. 2: The width of the switching current probability dis- 
tribution as a function of magnetic field with (solid line) and 
without (dashed line) taking into account zero-point oscilla- 
tions of CMs. Here we have taken typical parameters for the 
junction as I c o = WOfiA, and tiujpo/Ejo = 10 -3 , the junction 
size L = Aj. 



to the resistive state of the junction. In the quantum 
regime the process of escape is determined by the tunnel- 
ing of phase x an d the switching rate is given by [3, 0, S| 

48(v'2E J (ff)E /0 ) 1 / 2 / I a0 2 \ 5/4 

r(I,Q n )oce 5fi "p° v '^ H) ^™ anQ ") . 

(12) 

In experiment, the measured characteristics is usually the 
time-dependent probability P(I, t) of finding the junction 
in the zero-voltage state. Since, Q n are the random vari- 
ables, this probability can be written as 



P(I,t)= J p{Q n )exp[-tr(J,Q n )]dQ n 



(13) 



Here, p(Q n ) is the quantum-mechanical probability dis- 
tribution of CMs in the junction. 

If we neglect the excitation of CMs, Po(t) — 
exp(— t/ro) displays a standard exponential decay with 
time. Here, To = l/r(/,0) is the dc bias dependent life- 
time of the zero- voltage state, see Eq. (fT2"j). This behav- 
ior is shown in Fig. 3 by the dashed line. The deviation of 
P(t) from the exponential dependence allows to charac- 
terize the quantum-mechanical properties of CMs. The 
situation becomes especially interesting for nonequilib- 
rium and nonclassical states of the cavity modes. These 
are the zero-point oscillations, the chaotic, coherent and 
squeezed states, and various entangled states well known 
in quantum optics 0, Q • 

As a particular example, first, we consider zero-point 
oscillations in CMs. In this case the p{Q n ) is given as 
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Substituting Q14p in Eq. (|13[) and calculating all integrals 
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over Q n , we obtain for $ f» $0 
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where the magnetic field dependent parameter 71 — 

(t^t) 172 ( 2(S ) _1/4 - Fot sma11 ti mes t < H T o the 
influence of the equilibrium CMs quantum noise is small, 
and the dependence P(t) ~ exp(— t/ro) restores. Notice 
here, that for small external magnetic field and not very 
large number of excited CMs the exponential dependence 
of P(I, t) is always valid. 

Various nonequilibrium quantum-mechanical states of 
CMs can be induced by applying external microwaves. 
The interesting case is the coherent state of a single CM 
excited at a frequency w = u>\ . In this case all values of 
Q n except for Qi are small, and p(Qi) is given by 



PcohiQi) - exp 



(Qi-v) 2 



where 77 oc \W is determined by the power W of external 
microwave radiation. For a Josephson junction subject 
to an external magnetic field, we obtain that the exci- 
tation of the coherent state of a single CM leads to the 
dependence P(t) in the form 
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t > 72-Tcoh. (17) 



The parameter 72 



^"ejoS 5 ~ 1/2 de P ends on both 
magnetic field and the power of microwave radiation. 

Here r coh = 1/T(I,Q\ = i]), see Eq. ([T2]). Similarly 

to the case of a quantum noise in CMs, at short times 

t < 72T co h the P(I, t) decays exponentially with time, i.e. 

P (t) ~ exp(-t/r coh ). 

Another interesting quantum-mechanical state of CMs 

is the chaotic state of a single CM induced by microwave 

radiation 0, Q. In this case the power W of external 

microwave radiation determines only the mean photon 

number m, i.e. fh oc W. As the mean photon number 

m is relatively large the probability distribution p{Q n ) 

takes the form 
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Therefore, a strongly excited chaotic state of a single CM 
should show the same time-dependent probability P(t) 
as in the zero-point oscillations case (|15p , for which 71 is 
replaced by the microwave power dependent parameter 
73 = Ji/fh. 

The probability P(t) of finding the junction in the zero- 
voltage state is shown in Fig. 3 for various states of CMs. 
We note that more complicated states as squeezed states 
or Fock states (e.g., N photons in one mode n = 1) can 
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FIG. 3: Time-dependent probability P(t) of finding the junc- 
tion in the zero-voltage state: without taking into account 
CMs (dotted line) and with taking into account quantum 
fluctuations induced by CMs (solid line) and coherent state of 
CMs (dashed line). The parameters are chosen as 5 — 5 1CP 3 , 

and the 



(16) Ic(H)/I c (0) = 0.2, hw p0 /Ejo = KT 4 , 7] = 1.5 10 



junction size L 

TO = T coh 



0.27Aj. For simplicity, we have taken 



be prepared by using pulsed technique and intrinsic non- 
linearity of cavity modes • The entanglement of Fock 
states will manifest itself by oscillations of P(t) depen- 
dence. For realistic values of junction parameters we ob- 
tain 71 ~ 72 = 7, and the deviations of P(t) from the 
exponential decay should be detectable experimentally. 

In conclusion, we have shown that the excitation of 
cavity modes in distributed Josephson junction or paral- 
lel arrays of junctions manifests itself by either enhance- 
ment or suppression of the escape rate from the supercon- 
ducting state, depending on applied magnetic field. This 
effect is due to a renormalization of the potential barrier 
for the escape which, in turn, depends on the quantum 
state of the cavity mode. The important characteristics 
of the cavity mode quantum electrodynamics, namely the 
probability distribution of the CMs p(Q) can be detected 
experimentally by measuring the temporal decay of the 
switching probability P(t) given by Eqs. (fT2|) and (fl3|) . 
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